Assessing HCH isomer uptake in Alnus glutinosa: implications for phytoremediation and microbial response

Although the pesticide hexachlorocyclohexane (HCH) and its isomers have long been banned, their presence in the environment is still reported worldwide. In this study, we investigated the bioaccumulation potential of α, β, and δ hexachlorocyclohexane (HCH) isomers in black alder saplings (Alnus glutinosa) to assess their environmental impact. Each isomer, at a concentration of 50 mg/kg, was individually mixed with soil, and triplicate setups, including a control without HCH, were monitored for three months with access to water. Gas chromatography–mass spectrometry revealed the highest concentrations of HCH isomers in roots, decreasing towards branches and leaves, with δ-HCH exhibiting the highest uptake (roots—14.7 µg/g, trunk—7.2 µg/g, branches—1.53 µg/g, leaves—1.88 µg/g). Interestingly, α-HCH was detected in high concentrations in β-HCH polluted soil. Phytohormone analysis indicated altered cytokinin, jasmonate, abscisate, and gibberellin levels in A. glutinosa in response to HCH contamination. In addition, amplicon 16S rRNA sequencing was used to study the rhizosphere and soil microbial community. While rhizosphere microbial populations were generally similar in all HCH isomer samples, Pseudomonas spp. decreased across all HCH-amended samples, and Tomentella dominated in β-HCH and control rhizosphere samples but was lowest in δ-HCH samples.


Experimental design and sample processing
For this study, experimental substrates were prepared by mixing clean soil and standard α-HCH, β-HCH and δ-HCH isomers to achieve a 50 mg/kg dry weight mixture of each, with a control soil sample prepared in the same way but without HCH.Next, two-year-old A. glutinosa saplings were planted in triplicate for all experimental and control variants and then placed outdoors in an unshaded site with instant access to water.The saplings purchased at the regular forestry nursery (Lesoškolky, Vestecká 999, 250 01 Brandýs nad Labem-Stará Boleslav, Czech Republic).This nursery is a licensed source of plant reproductive material.No further permissions are needed under Czech law to purchase saplings.The nursery declared the purchased saplings were A.glutinosa.We confirm that the experimental work performed with the plant material in this study complies with relevant institutional, national, and international guidelines and legislation.
After three months, 10.0 ± 0.5 g of soil was sampled from each pot for analysis of α-HCH, β-HCH and δ-HCH isomers, with an equal amount of soil used for dry mass determination (see below).The A. glutinosa saplings were then removed from the pots and divided into four sections, i.e. root, trunk, branches and leaves.Each part of the sapling was then weighed, after which it was milled in liquid nitrogen and the resultant fine powder used for HCH analysis, dry weight determination (both using 5.0 ± 0.5 g of biomass) and phytohormonal analysis (see below).
Samples for HCH extraction were agitated in a horizontal shaker at 200 rpm for 24 h in 10 ml of a 1:1 (v/v) acetone:hexane mixture for soil, and 5 mL of 1:1 (v/v) acetone:hexane for sapling biomass.The supernatant was then collected and anhydrous sodium sulphate added, after which the mixture was vortexed.One ml of extract was then taken and enriched with 10 μl of internal deuterated γ-HCH standard to cover ionisation instability prior to analysis (see below).Samples for dry mass determination were dried for three hours at 105 °C and then reweighed.
Rhizosphere soil, obtained by shaking the roots, and any remaining soil stuck to the roots was placed in a 50-ml test tube containing 35 ml of sterile phosphate buffer and then shaken for 5 min.The resultant solution was then centrifuged for 10 min at 4000×g at 4 °C to obtain biomass for DNA extraction (see below).

Chemical analysis
Prior to mixing the experimental substrates, the purity of the standard isomer solutions was assessed using a RSH/Trace 1310/TSQ8000 Gas chromatography-tandem mass spectrometry (GC-MS/MS) assembly with a DB-5 ms column (PAL, Switzerland; ThermoFisher Scientific, USA).The same assembly was also used to assess the concentrations of HCH and ClB isomers in the various experimental samples.In each case, the limit of quantification (LOQ) was set at < 0.01 µg/g (dw) for all analytes.The standard deviation is up to 20% due to the complexity of the biological matrix.
Phytohormone analysis was performed using an Acquity® I-Class ultra-high-performance liquid chromatograph (UHPLC; Waters, Milford CT, USA) coupled with a Xevo TQ-XS MS/MS assembly (Waters, UK), using the isotope dilution method described by Šimura et al. 54 .Values under the limit of detection (LOD) were replaced by the values representing 95% of the detection limit for analysis.

DNA extraction and real-time quantitative PCR
DNA extraction of soil and rhizosphere samples was undertaken in duplicate using the DNeasy power Soil KIT (Qiagen, Netherlands).DNA yield and quality were assessed using a Qubit fluorometer (Thermo Fisher Scientific, USA) and agarose gel electrophoresis.

Amplicon 16S rRNA sequencing
The V4 region of the bacterial 16S rDNA gene was amplified using the primers 530F (5′TGC CAG CMGC-NGCGG3′) and 802R (5′TACNVGGG TAT CTA ATC C3′), while for fungal abundance, the ITS2 region was amplified using the primers ITS3-F (5′ GCA TCG ATG AAG AAC GCA GC3′) and ITS4-R (5′ TCC TCC GCT TAT TGA TAT GC3′) at a final volume of 50 µL 59,60 .The amplicons were cleaned using the Agencourt Ampure XP system (Beckman Coulter, USA) and the Ion Torrent platform (Thermo Fisher Scientific, USA) was used for sequencing analysis.Barcoded sequencing adapters were ligated to the PCR products using the Ion Xpress Plus gDNA fragment library kit with Ion Xpress barcode adapters (Thermo Fisher Scientific, USA) and the samples analysed using the Ion PGM Hi-Q Sequencing Kit using an Ion 314 Chip (Thermo Fisher Scientific, USA).
The raw Ion Torrent reads were first processed using QIIME 2 v.2021.8software 61 , after which the raw sequence data were demultiplexed and quality filtered using the q2-demux plugin followed by denoising with DADA2 62 , with reads below 250 bp being removed.Taxonomy was assigned to each amplicon sequence variant (ASV) using the q2-feature-classifier 63 classify-sklearn naive Bayes taxonomy classifier against the Silva 138 database 64 , after which Mitochondria and Chloroplasts were removed.Accuracy of classification was evaluated against an artificial MOCK community sample.The QIIME 2 outputs were then processed using the phyloseq R package 65 .

Statistical analysis
The effect of each HCH isomer on sapling growth parameters was evaluated using one-way ANOVA with posthoc Tukey tests, using the Origin software package v.2019b (OriginLab Corporation, USA).Prior to analysis, all data were subjected to Levene's test to assess homogeneity of variance.All statistical analyses were performed at a significance level of α > 0.05.Data for phytohormonal analysis were first log-transformed and then subjected to principal components analysis (PCA) to describe the general data structure.The non-parametric Mann-Whitney test was then used to calculate significant differences in phytohormonal content between the control and HCH-polluted saplings.

Sapling root biomass
The root biomass of saplings treated with α-HCH contained 7991.2 ng/g of the isomer, and a total concentration of 898.5 ng/g ClB, both metabolites matching those found in the soil.A similar decrease in metabolite concentrations was also observed for the β-HCH isomers, with root concentrations reaching 1134.63 ng/g and total metabolite levels just 80.43 ng/g (primarily hexaClB, pentaClB and 1,2,4-triClB).

Above-ground sapling parts
In general, concentrations of HCH were low in trunks, branches and leaves.Nevertheless, saplings exposed to α-HCH had concentrations reaching 2677.16ng/g in trunks, 9839.33 ng/g in branches and 1506.30ng/g in leaves, indicating that, once taken up by the roots, the isomer is subsequently transported throughout the sapling (Fig. 1).Interestingly, highest concentrations were determined in branches, where, presumably, no further processing of the original substance took place as no metabolites were detected.Subsequently, the isomer is transported to the leaves where further metabolism or volatilisation does occur.It should be noted, however, that there is a possibility that hitherto unknown metabolites are formed in all sapling parts.Overall, less HCH and metabolites were recorded in sapling parts treated with the β-HCH isomer, with levels reaching just 273.36 ng/g in trunks, 560.60 ng/g in branches and 701.30ng/g in leaves.Likewise, total metabolite concentrations were negligible, with just 1,2-and 1,4-diClB and 1,2,4-triClB detected in the trunks and hexaClB in the branches.
Results for the δ-HCH isomer were like those for the α-HCH isomer, with trunks containing 1108.63 ng/g, branches 5469 ng/g and leaves 1916.30ng/g.Once again, there were higher levels in the branches than in other aboveground parts, though metabolites were only detected in trunks and not in branches or leaves.
Phytohormones.Overall, 72 plant hormones and metabolites in eight classes were identified in tissues growing in contaminated and control soils (Fig. 2).The roots, for example, were characterised by highest concentrations of 2-methylthio CK and indole-3-acetamid.In comparison, highest concentrations of IAA (bioactive AUX) and the majority of CK were detected in sapling trunks.Branches and leaves typically had high concentrations of BR, accompanied by increased levels of ABA, JA and SA.CK in all sapling organs was downregulated in trees growing on polluted soil, whatever the isomer (Fig. 3).The saplings also contained metabolites of GA, with a significant increase in the precursor GA53 and the catabolite GA29 in stressed trees, along with upregulated JA and downregulated abscisates (Fig. 3a).S1.The ellipse expresses the 95% limit of multivariate t-distribution (Hotelling's T-squared distribution).The similarities of the samples are given by the concentrations of plant hormones depicted in the loading plot (b).The plant hormones co-localising with the samples on the same side of the scatter plot are abundant in these samples.Conversely, the plant hormones on the opposite side of the scatter plot are low in their concentration.The plant hormones are coloured by the corresponding phytohormonal group.

Bacteria
Soil and rhizosphere samples produced 530 common bacterial taxa.Of these, 93 taxa were found in the rhizosphere but not in soil, and 126 in soil samples only (Fig. 1 supplementary).Soil bacterial community structure was similar in all samples treated with HCH.
The dominant bacterial taxa in soil samples were, in order, Pseudomonas, the unclassified Chloroflexi group KD4-96, Comamonadaceae, Nocardioides and Gemmatimonadaceae.While Pseudomonas was dominant, with abundance levels highest in control samples, levels decreased noticeably in samples treated with HCH (Fig. 4).
Bacterial abundance was similar in all rhizosphere population samples, though with some exceptions (Fig. 4).For example, rhizosphere Pseudomonas levels were significantly lower than those in soil for all HCH treatments, with lowest abundance for the δ-HCH isomer, despite being detected at highest levels in the root biomass (Fig. 1).In contrast, Pseudolabrys, Gaiellales (Gaiella), Acidobacteriales, Haliangium and Bauldia were the dominant contributors to biomass in the rhizosphere.

Fungi
Soil and rhizosphere samples contained 138 common fungal taxa, with 33 taxa found in the rhizosphere but not in soil and 86 found in soil only (Fig. 1 supplementary).The abundance of Penicillium and Coleophoma was higher in the soil than in the rhizosphere, while levels of Mortierella were higher in the rhizosphere (Fig. 5).While Tomentella was generally dominant in rhizosphere samples, samples treated with δ-HCH had the lowest concentrations and β-HCH samples and the control the highest.

Functional genes involved in HCH biodegradation
The lin genes were found in all samples, with higher relative quantities in soil than the rhizosphere (Table 1).Interestingly, the linA gene was not detected in most rhizosphere samples but was present at low quantities in all soil samples.The linB gene was detected at high levels in soil samples, especially in samples treated with β-HCH.Overall, highest levels of the linB-RT gene were found in soil samples treated with the β-HCH isomer, while very little was recorded in control samples and none at all in soil treated with α-HCH, suggesting that the β-HCH isomer biodegrades in soil 44 .Finally, linD genes were detected at intermediate relative quantities in all rhizosphere samples, but were rarely detected in soil samples.

Discussion
This study provided several important findings regarding the impact of HCH pollution on A. glutinosa saplings and its associated soil microbiota, with implications for understanding the effects of persistent organic pollutants (POPs) on plant growth, microbial communities,phytohormone regulation, and applicability of HCH concentration in tree biomass for elucidation of HCH subsurface contamination (HCH phytoscreening).

HCH removal efficiency
HCH removal efficiency varied with type of HCH isomer, with β-HCH exhibiting highest removal efficiency (90.26%), followed by α-HCH (64.85%) and δ-HCH (57.08%).Surprisingly, a high amount of isomer β-HCH was removed in this experiment, contrary to its persistent properties.In contrast, the removal of 64.58% of α-HCH is consistent with the findings of the study by Liu et al. 66 , where 47.3% was removed under hydroponic conditions in the presence of wheat.Chen et al. 67 51 may be due to young age of the A. glutinosa seedling in their experiment.The differences in presented results can be attributed to differences in the persistence and physical and chemical properties of each HCH isomer, such as solubility and volatility.As the δ-HCH group is the most lipophilic (log K ow = 4.14), soil composition may have negatively affected the availability of this isomer group to saplings and microbiota 69 .Simultaneously, this significant decline may be caused by the environmental influences.These findings were demonstrated in a study by Liu et al. 70 , where, when comparing the initial β-HCH concentration with the final concentration in the unplanted variant, this concentration was reduced by 50%.It is important to acknowledge, however, that the seedlings in this study were exposed to freshly dissolved HCH isomers, while at actual contaminated sites, these isomers may undergo chemical transformation due to interactions with environmental factors such as soil, water and air.

Metabolite levels in HCH polluted soils
Lowest ClB metabolite levels were recorded in β-HCH polluted soils.It is important to note that transportation of HCH does not always follow the 'water diffusion pattern' through the environment as its octanol-water partition coefficient is close to 4 (log K ow ) 69 .In contrast, soils containing δ-HCH had 16.7 ng/g of ClB metabolites, primarily consisting of the pentaClB and hexaClB isomers.In comparison, α-HCH is more susceptible to metabolic transformation in soils, most likely through bacterial activity in cooperation with plants.Indeed, numerous studies have reported high degradation rates of α-HCH through a range of microorganisms [71][72][73] .
Based on the available literature, it is evident that ClB can be metabolized by microorganisms and possibly by plants.Trapp et al. 74 demonstrated the presence of polar metabolites of ClB in plants, indicating their potential metabolism by plants.Furthermore, Middeldorp et al. 75 and Kurt & Spain 76 provided evidence of ClB biotransformation under specific conditions, suggesting microbial metabolism.The studies by Wang et al. 77 and Langenhoff et al. 78 also indicated the biodegradation potential of ClB by microbial consortia under different conditions.Moreover, Ramanand et al. 79 and Nishino et al. 80 provided insights into the reductive dehalogenation and degradation of ClB by bacteria, further supporting the microbial metabolism of these compounds.

HCH uptake and transport in saplings
The highest concentration of HCH isomers was recorded in the root part of plants.This is attributed to the proximity of contaminated soil and to the lipophilicity of the individual isomers and the root structures.The lipophilic compounds are able to enter the root tissue, where in most cases they remain bound to the lipophilic layer of the root epidermis 81 .The uptake of HCH by plants has been recorded under laboratory 70 and field conditions 16 .The root biomass of saplings treated with α-HCH contained 7991.2 ng/g α-HCH and a total ClB concentration of 898.50 ng/g, mirroring soil metabolite levels, and exhibited a rapid decrease in metabolite concentrations, possibly due to factors such as low uptake, rapid degradation, transportation or transformation within the roots.A similar decline in metabolite concentrations was also evident for the β-HCH isomers.In terms of above-ground sapling components, HCH concentrations were generally low in trunks, branches, and leaves.Saplings exposed to α-HCH, however, exhibited elevated concentrations in leaves, suggesting that once absorbed by the roots, the isomer is subsequently transported throughout the sapling.Intriguingly, substantial concentrations were identified in branches, where it is presumed that minimal further processing of the original substance occurs as no metabolites were detected.Conversely, the isomer is then transported to the leaves, where additional metabolism or volatilization may take place.This transport is supported by the vapor pressure of the isomer, which may explain the transport from the root to the aerial parts of plants 82 .It is essential to acknowledge the possibility of the formation of yet undiscovered metabolites within all sapling parts.In contrast, sapling parts treated with the β-HCH isomer displayed lower HCH and metabolite levels, while the results for δ-HCH closely resembled those observed for α-HCH.Similarly, higher concentrations were identified in the branches compared to other aboveground parts, though ClB metabolites were only detected in the trunks and not in branches or leaves.

Phytohormonal responses
The greatest difference in phytohormone concentration levels was observed between the various sapling organs, most likely due to differences in organ function and developmental patterns over time and space.For example, the roots showed notable levels of 2-methylthio CK and indole-3-acetamid, which serve as indicators of bacterial activity in the rhizosphere 83,84 .The highest concentrations of IAA (bioactive AUX), and the majority of CK, were detected in sapling trunks, presumably as homeostasis between trunk AUX and CK regulates various physiological functions, including apical dominance 85 .The phytohormonal profile of branches and leaves displayed enrichment of BR, typically associated with the meristems of actively growing young organs 86 .Additionally, elevated levels of ABA, JA and SA in the branches and leaves indicated a potential response to abiotic stress conditions 87 .Across all sapling organs, CK levels were generally downregulated in trees cultivated in contaminated soil, irrespective of the HCH isomer applied.This downregulation may be linked to early senescence as a reaction to stress induced by phytoremediation efforts 88 .Notably, isopentenyladenine and isopentenyladenosine in the sapling trunks showed exceptions to this trend, likely associated with CK transport mechanisms 89,90 .Furthermore, our study identified metabolites of GA in the A. glutinosa samples, which have previously been described as components of the 13-hydroxy biosynthetic GA pathway found in buds of mountain alder Alnus incana subsp.tenuifolia 91,92 .In this study, stressed trees exhibited a significant increase in precursor GA53 and the catabolite GA29, alongside upregulated JA and downregulated ABA.The notable concentration of JA suggests that signalling pathways involved in abiotic stress responses had been activated.It is noteworthy that components of the JA signalling pathway have previously been reported as interacting antagonistically with ABA in Arabidopsis (rockcress, Brassicaceae), which could explain the observed lower concentration of ABA metabolites in our study

Microbial communities in soil and rhizosphere
Bacterial community analysis revealed a comprehensive array of 530 common bacterial taxa in both soil and rhizosphere samples.Furthermore, our findings revealed the prevalence of Nocardioides, Chloroflexi and Gemmatimonadaceae, which have also been reported as abundant taxa in soils contaminated with organochlorine pesticides (OCPs) 96 .Interestingly, the unclassified Chloroflexi group KD4-96 has frequently been detected in soils contaminated with heavy metals 97 , with abundance correlating strongly with iron (Fe) and aluminium (Al) concentrations 98 .In contrast, Pseudolabrys, Gaiellales (Gaiella), Acidobacteriales, Haliangium and Bauldia were all dominant contributors to biomass in the rhizosphere, with Pseudolabrys, order Rhizobiales, being nitrogenfixing bacteria that establish a symbiotic relationship with plant roots, similar to Bauldia.Additionally, Gaiella plays a crucial role as an organic matter decomposer and actively participates in carbon cycling 99 .Of the fungi, Tomentella was generally dominant in rhizosphere samples, while soil samples treated with δ-HCH had very low concentrations and soils with β-HCH and the control group having highest levels.Tedersoo et al. identified 40 species of putatively ectomycorrhizal fungi from seven sites dominated by A. glutinosa, with the ectomycorrhizal fungi Tomentella aff.Sublilacina most prevalent under saline stress conditions 100 .Tomentella contain melanins, which may act as a boundary between fungal cells and their environment, protecting them against physical, chemical and/or biological stressors 101 .Our findings revealed that five taxonomic bins: Ferribacterium, Methylotenera, Fluvicola, env.OPS_17 and Exiguobacterium exhibited statistically significant differences exclusively within soil samples between experimental treatments and control counterparts (Fig. 7 supplementary and Table S2).Conversely, in other sample types, the observed distinctions did not attain statistical significance.

Functional genes related to HCH biodegradation
Several functional genes related to HCH biodegradation were detected, including linA, linB and linD.The lin genes were found in all samples, with higher quantities in soil than in the rhizosphere.Presence of the linB gene was most pronounced in soil samples, particularly in those exposed to β-HCH.Lal et al. (2010) also reported linB occurring at higher levels than linA in most soil samples from a highly contaminated HCH dump site 44 .The presence of linB genes in soils with β-HCH may indicate the potential for efficient biodegradation of HCH isomers, aligning with the observed relationship between linB genes and β-HCH removal efficiency.Sharma  et al. (2006) demonstrated that Haloalkane Dehalogenase LinB is responsible for β-and δ-HCH transformation in Sphingobium indicum B90A, highlighting the enzymatic basis for the degradation of HCH isomers, including β-HCH.The key function of linB is the breakdown of the most resistant HCH isomer, β-HCH 102 .Furthermore, the similarity of linB to other enzymes (as opposed to linA, which is unique) may result in the selection of more linB genes from unidentified bacteria present in the samples 44 .These findings suggest that soil microbial communities possess the genetic potential for HCH biodegradation, with linB a prominent gene involved in β-HCH degradation.The presence of these genes in the rhizosphere indicates the potential for plant-microbe interactions in HCH remediation.

Study limitations
We recognise that the absence of an unplanted control group in this study may be a limitation.While we endeavoured to investigate the bioaccumulation potential of HCH isomers in A. glutinosa saplings and their associated microbial communities comprehensively, the inclusion of an unplanted control group would have provided valuable insights into the background levels of HCH decontamination via exterior physical conditions.Despite this limitation, our study offers significant insights into HCH isomer bioaccumulation patterns in relation to A. glutinosa saplings and their associated responses, contributing to the broader understanding of HCH persistence in the environment and its ecological implications.

Discussion of hypotheses
Our first hypothesis, which suggested that the δ-HCH isomer would show lowest bioaccumulation in A. glutinosa saplings, followed by α-HCH and β-HCH, was not supported by the findings.The δ-HCH isomer did show, in contrast, highest bioaccumulation potential, followed by α-HCH and β-HCH.These differences in bioaccumulation can be attributed to variations in chemical properties and affinities of individual HCH isomers for plant uptake.
Our second hypothesis proposed that highest HCH concentrations would be found in the root system of A. glutinosa saplings, with decreasing concentrations in the trunk, branches and leaves, reflecting the plant's uptake and translocation mechanisms.This hypothesis was largely confirmed by our results, which showed that HCH concentrations were indeed highest in the roots, with sequentially decreasing levels in above-ground plant parts.This distribution pattern aligns with the plant's mechanisms for HCH uptake and translocation.
Our final hypothesis suggested that the introduction of HCH isomers would lead to changes in the composition and abundance of microbial communities, with specific bacterial species responding differently to different HCH isomers, potentially indicating their ability to degrade or interact with the contaminants.The observed microbial community dynamics in response to HCH isomers were either subtle or not statistically significant, suggesting that specific populations may not have exhibited pronounced variations in their response to different HCH isomers.This outcome implies that factors beyond HCH isomer type, such as soil conditions or other environmental variables, may have played a more influential role in shaping microbial communities during the study.Further research including a more comprehensive exploration of these factors will be necessary to www.nature.com/scientificreports/improve our understanding of the intricate interactions between microbial communities and HCH isomers in contaminated environments.

Conclusion
A.glutinosa saplings showed clear uptake of all HCH isomers, with highest quantities detected in the roots and lowest in the leaves.Owing to its physicochemical properties, the δ-HCH isomer was the most persistent in soils and the most strongly bound to A. glutinosa roots, with α-HCH the second major isomer recorded in soils treated with β-HCH and δ-HCH.All metabolites were present at highest concentrations in the α-HCH isomer treatment, suggesting that α-HCH may be subject to degradation over the longer term.All HCH isomers were found at highest proportions in the soil, with relatively little found in root biomass, suggesting that degradation of HCH isomers by soil bacteria occurs mainly through the upstream pathway.However, high levels of the linD gene in rhizosphere bacteria, also confirmed HCH degradation via downstream pathways.Overall, there was no significant difference in the abundances of bacterial and fungal consortia between treated and control samples.Similarly, there were no significant differences between soil and rhizosphere microorganisms.Phytohormone analysis indicated that A. glutinosa reacts to HCH contamination through changes in the stress hormones CK, JA, abscisate and GA.
Our findings have broad practical implications for a range of fields.Firstly, they offer valuable insights for environmental remediation efforts, aiding in the development of more effective and sustainable phytoremediation strategies for HCH-contaminated soils.Secondly, this research contributes to the assessment of soil and ecosystem health in areas affected by POPs, thereby informing ecosystem restoration and conservation strategies.Moreover, these insights could drive advancements in bioremediation technologies, offering naturebased solutions for contaminated land restoration.While our results provide preliminary information on the phytoremediation potential of A. glutinosa trees, further studies will be needed at actual contaminated sites to fully understand their ability to remediate contaminated soils. 1506

Figure 2 .
Figure 2. Variability of A. glutinosa saplings in plant hormone concentrations calculated by principal component analysis (PCA).In the score plot (a), each data point represents an individual plant sample coordinated by the principal components PC1 and PC2, both expressing total variability 45.3%.The similar samples, each coloured by a plant organ, are clustered close to each other.The numbers nearby the data points express the sample identifiers, whose meanings are explained in TableS1.The ellipse expresses the 95% limit of multivariate t-distribution (Hotelling's T-squared distribution).The similarities of the samples are given by the concentrations of plant hormones depicted in the loading plot (b).The plant hormones co-localising with the samples on the same side of the scatter plot are abundant in these samples.Conversely, the plant hormones on the opposite side of the scatter plot are low in their concentration.The plant hormones are coloured by the corresponding phytohormonal group.

Figure 3 .
Figure 3.Effect of HCH on phytohormone concentration levels in the different parts of A. glutinosa; (a) branches, (b) trunk, (c) leaves and (d) roots.The blue circles in each volcano plot represent the plant hormones distributed according to the mean rank difference of the concentration between control (C) and HCH-treated groups (horizontal axis) and the negative logarithm of p-value calculated by Mann-Whitney test (vertical axis).The horizontal dashed line delimits the significance level of α = 0.05.The blue boxes inside each volcano plot shows the data distribution of the most significantly different plant hormone concentrations (with the highest absolute value of both mean rank difference and negative logarithm of p-value) in each individual sapling.The asterisks indicate the statistical significance at level of α = 0.05 (p-value 0 '****'0.0001'***'0.001,'**'0.01).